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ABSTRACT 

We present a survey toward 16 Low Luminosity Objects (LLOs with an internal luminosity, Lint) 
lower than 0.2 Lq) with N2H+ (1-0), N2H+ (3-2), N2D+ (3-2), HCO+ (3-2) and HCN (3-2) using 
the Arizona Radio Observatory Kitt Peak 12m Telescope and Submillimeter Telescope. Our goal is to 
probe the nature of these faint protostars which are believed to be either very low mass or extremely 
young protostars. We find that the N2D+/N2H+ column density ratios of LLOs are similar to those of 
typical starless cores and Class 0 objects. The N2D+/N2H+ column density ratios are relatively high 
(> 0.05) for LLOs with kinetic temperatures less than 10 K in our sample. The distribution of N2H+ 
(1-0) line widths spreads between that of starless cores and young Class 0 objects. If we use the line 
width as a dynamic evolutionary indicator, LLOs are likely young Class 0 protostellar sources. We 
further use the optically thick tracers, HCO+ (3-2) and HCN (3-2), to probe the infall signatures of 
our targets. We derive the asymmetry parameters from both lines and estimate the infall velocities by 
fitting the HCO’*' (3-2) spectra with two-layer models. As a result, we identify eight infall candidates 
based on the infall velocities and seven candidates have infall signatures supported by asymmetry 
parameters from at least one of HCO+ (3-2) and HCN (3-2). 

Subject headings: star: brown dwarfs - stars: low-mass - stars: protostars - ISM: kinematics and 
dynamics - ISM: molecules - radio lines: ISM 


1. INTRODUCTION 

Very Low Luminosity Objects (VeLLOs) are embed- 
ded protostars with intern al luminosity Lint < 0.1 Lq 
(|Di Francesco et al. I 2007tl. The first VeLLO, L1014- 
IRS, was identihed bv lYoung et al.l (I2004D toward a dense 
core previously thought to be starless, L1014. The 
SED fitting results suggest that L1014-IRS has a very 
low internal luminosit y ^ 0.09 Lit, contributed by the 
centr al star and disk (jYoung et al.l 12004 iHuard et al. I 
I2OO6II . Three possibilities are proposed to exp lained 
VeLLOs’ low luminosities (|Dunham et al. II2014D : they 
can be 1) very low mass protostars, 2) extremely young 
protostars, or 3) protostars in a quiescent phas e of an 
episodic accretion process. lAndre et al. I (1201211 identi¬ 
fied a “pre-brown dwarf\ Oph-Bll, and suggested that 
brown dwarfs can form in iso lated cores as low mass stars 
(|Padoan &: Nordlund Il2004l . By comparing with evolu¬ 
tionary track from models, the low luminosities imply 
that VeLLOs are substellar sources (current mass < 0.08 
Mq) embedded in molecular cores (L1014-IRS: Young et 
al. 2004; Huard e t al. 2006, L1521F: Bourke et al. 2006). 
iLee et aL~l (|2013H suggested that the mass of L328 is at 
most 0.05 Mq based on the estimate of mass loss rate 
from outflows. Therefore, VeLLOs could form low mass 
stars or even brown dwarfs depending on their future 
accretion. VeLLOs with small mass accretion reservoir 
(parent cores) are very likely brown dwarfs still in the em¬ 
bedded phase, i.e. “proto brown dwarfs’’ (L328: Lee et 
al. 2009, 2013, L1148: Kauffmann et al. 2011, ICE 348- 
SMM2E: Palau et al. 2014). Since VeLLOs are usually 
isolated, they could be brown dwarfs forming through 
cloud fragmentation followed by gravitational collapse. 
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On the other hand, VeLLOs’ substellar masses also 
indicate that VeLLOs can be early Class 0 protostars 
or even younger than Class 0 protostars (IRAM 04191: 
Andre et al. 1999; Belloche et al. 2002; Dunham et al. 
2006, L1521F: Bourke et al. 2006; Takahashi et al. 2013, 
Cha-MMSl: Belloche et al. 2006; Tsitali et al. 2013; 
Vaisala et al. 2014). Co nsidering seve r al age indicators 
such as Tboi and Lboi, lAndre et al. I (|1999f) suggested 
that IRAM 04191 is a young Class 0 object with an age 
of < 3-5 X 10^ yr. iBelloche et al. I (j2002ll later con¬ 
firmed this by considering the inner free fall region ex- 
panding at the sound s peed. By studying the outflows, 
iTakahashi et al. 1 (l2013f) found that LI52IF is probably 
a pro tostar in the ear l iest ev olutionary stage (age < 10^ 
yr). IBelloche et al. I (|2006D found a high level of deu¬ 
terium fraction and no sign of large scale outflows for the 
VeLLO Cha-MMSl, suggesting that it has not driven any 
outflow yet and could be younger t han IRAM 04191 and 
L1521F. The kinematical study in iTsitali et al. I (j2013f ) 
further supports that Cha-MMSl is either a you ng Class 
0 sou rce or even a first-hydrostatic core (FHSC) (iLarson I 
Il969f) . The FHSC phase represents the short-lived phase 
during which a hydrostatic compact object forms at the 
center of a collapsing dense c ore but has no t yet experi¬ 
enced the dissociation of H2. (jLarson Ill969| l. The disso¬ 
ciation of H2 at the end of the FHSC phase triggers the 
second collapse that leads to the formation of the proto¬ 
star. Several FHSC candidates have been identified but 
none has been firmly confirmed yet (Cha-MMSl: Bel¬ 
loche et al. 2006, 2011a, L1448-IRS2E: Chen et al. 2010, 
Per-Bolo 58: Enoch et al. 2010; Dunham et al. 2011, 
L1451-mm: Pineda et al. 2011, CB17-MMS: Chen et al. 
2012, Bl-bS, and Bl-bN: Pezzuto et al. 2012, Murillo & 
Lai 2013). 
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The third possibility - a protostar in a quiescent ac¬ 
cretion phase - is possible because the mass accretion 
can enhance the luminosity of protostars, i.t. the ac¬ 
cretion luminosity Lace (L673-7: Dunham et al. 2010a, 
CB 130: Kim et al. 2011 ) . Episodic accretion mod- 
els (|Dunham et al. l[2010bl : iDiinham h Vorobvov1l2012D 

have been proposed to interpret the low luminosities of 
VeLLOs and could be a possible s olution of the luminos¬ 
ity problem (|Kenvon et al. I[l990f) . The accretion lumi¬ 
nosities estimated based on observed CO outflows further 
support that th e accretion process is episodic i n low lumi- 
nosity objects (|Dunham et al. 1 12006112010at ILee et al~1 
[Moll . Note that these three possibilities do not conflict 
with each other; VeLLOs could be a m ixture of these 
three possibilities (iDunham et al. Il2014ll . 

The most tho r ough VeLLO survey has been done by 
IDunham et al. I (j2008H . The Spitzer Legacy Project 
“From Molecular Cores to Planet Forming Disks” (c2d; 
Evans et al. 2003, 2009) mapped five nearby molec¬ 
ular clouds and 95 small dense cores, which provides 
an excellent opportu nity for searching for VeLLOs. 
IDunham et al. I (j2008ll found a tight relation between 
the 70 /im flux of protostars and their internal luminos¬ 
ity Lint = 3.3 X 10® Lq based on 1460 synthetic 
SEDs fro m the Monte Carlo dust radia tive transfer code 
RA DMC (iDullem ond fc Dominildl2004ll . Using this rela¬ 
tion, jOunha^^^OlloO^MMructed a list of embed¬ 
ded protostars with low Lin t and identified 15 V eLLOs 
from the c2d data. Lately, iHsieh fc Lai I (l2013| l devel¬ 
oped a new YSO identification method to probe faint 
YSOs and identified seven additional VeLLO candidates 
from the c 2 d survey clouds. 

In order to understand the nature of VeLLOs, we focus 
on probing their evolutionary states using chemical and 
dynamic evolutionary tracers. It is known that at the 
center of starless cores, where the temperature is about 
10 K and the number density is above 10^ cm“®, common 
molecular trace rs such as CO and CS freeze out onto dust 
grain surfaces ([Caselli et al. 1 119991 : iTafalla et al. I I2002L 
120041) . Other tracers are required to investigate the dy¬ 
namical and chemical properties in the cold dense region. 
The nitrogen-bearing species such as N 2 H+ and NH 3 sur¬ 
vive in the gas p hase at least for densities in the range 
of 10 ®— 10 ® cm~® (iCrapsi et al. ll2005bHJohnstone et al. I 

l2010HRoberts fc Millar 11200 71) . Thus, they are important 

tracers during the prestellar and protostellar phase. 

The N 2 D+/N 2 H+ abundance ratio can be an evolu¬ 
tionary ind icator in early evolu t ionary stages of star 
forma tion (iCraosi et al. I l2005bt lEmprechtinger et al. I 
120091) . The N 2 D+/N 2 H+ abundance ratio is affected 
by the exothermic reaction, H;^- 1 -HD ^ H 2 D+-I- H 2 , 
since H oD+ transfers the deuteron to N 2 and produces 
N 2 D+ ( Pore et al. 1 1200 J) . Because CO is a destroyer 
of H 2 D+ and (jCaselli fc Ceccarelli 1 12012}) . the de¬ 
pletion of CO increases the abundance of H 2 D+ and H)|’. 
The abundance increase of the latter molecule also speeds 
up the reaction to produce H 2 D+. As a result, the CO 
depletion increases the abundance of H 2 D+ as well as 
N 2 D+ and leads to a high N 2 D'^/N 2 H+ abunda nce ra¬ 
tio in a cold dense region. iCrapsi et al. I (l2005bD found 
that, in 31 low mass starless cores, the N 2 D+/N 2 H+ col¬ 
umn density ratio (N(N 2 D+)/N(N 2 H+) ratio) is corre¬ 
lated with several evolutionary indicators such as core 


density, line width and line asymmetry, suggesting that 
the N(N 2 D^)/N(N 2 H'*~) ratio increase s as a starless core 
evolves. lEmprechtinger et al. I (|2009D showed that the 
N(N 2 D+)/N(N 2 H+) ratio decreases as the core evolves 
while getting hotter and less dense after the central 
protostar formed. Thus, the N(N 2 D"*')/N(N 2 H+) ratio 
increases in the prestellar c ore phase and declines in 
the protostellar core p hase (jMiettinen fc Offner 1 120131 : 
iHuang fc Hirano1l2013D . 

Infall motions are also used as an evolutionary indi- 
cator since they are likely t o occur in the early stages 
([Emprechtinger et al. 1 120091 ) . Infall motions in proto¬ 
stars can be traced with observations including both op¬ 
tically thin and thick lines. For a protostar with an en¬ 
velope dominated by infall, the foreground gas and back¬ 
ground gas are red-shifted and blue shifted, respectively, 
with respect to the systemic velocity. If the foreground 
gas is colder than the central region, the red-shifted emis¬ 
sion of an optically thick line can be absorbed, resulting 
in a spectral line asym metry with a blue-shifted peak 
(iMardones et al. 1I1997D . 

We study the chemical and physical properties of 16 
Low Luminosity Obje cts (LLOs, Lint < 0.2 Lq) from 
IDunham et al. I (j2008f ) . We describe our sample and ob¬ 
servations in §2 and report the results in §3. The data 
analysis is presented in §4 and the properties of LLOs 
are discussed in §5. We summarize the results in §6. 

2. OBSERVATIONS 
2.1. Sample 

Using the c2d d ata ([Evans et al.l 120031 I2009D . 
IDunham et al. I (I2008D identified 50 embedded protostars 
with Lint < 1 Lq and 15 sources out of them are VeLLOs 
(Lint < 0.1 Lq). Hereafter we use the combination of 
the first letter of the three authors’ last names plus the 
source number in their paper as our source name {e.g. 
DOE 185). Our sample consists of 10 VeLLOs out of the 
15 sources. We excluded five sources: 1) DCE145 was 
identified as a background galaxy based on our CEHT Ks 
band (2.146 /im) outflow survey (Hsieh & Lai in prepara¬ 
tion), 2) DOE 161 and DOE 018 are unobservable using 
telescopes at the Arizona Radio Observatory (ARO), and 
3) DOE 024 (CB-130) and DOE 032 (L1148) have spec¬ 
tra contaminated by their reference positions that were 
by mistake set too close to the sources. In order to re¬ 
duce the bias from the small sample, we added 6 Low 
Luminosity Objects (LLOs, 0.1 Lq < Lint < 0.2 Lq) 
located in the two closest molecular cl ouds. Ophiuchus 
(125 p c) and Perseus (250 pc) from the IDunham et al. I 
(|2008l) sample. This sample includes 10 Class 0 sources 
and 6 Class I sources (bolometric temperature Tboi > 70 
K, Chen et al. 1995) as shown in Table 1. These sources 
show evidence that they are indeed embedded, allowing 
us to study the chemical and physical properties of their 
parent cores. 

2.2. Kitt Peak 12m Telescope observations 

We observed N 2 H+ (1-0) toward the 16 LLOs in 
November 2012 and February 2013 using the Kitt Peak 
12 meter Telescope (KP12m) at the Arizona Radio Ob¬ 
servatory (ARO). We used the position -switching mode 
point ed toward the infrared sources ([Dunham et al. I 
I 2 OO 8 D (Table 1). The integration time (including on and 
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TABLE 1 
LLO SAMPLE 


Source number 

Other name 

R.A. 

(J2000.0) 

Dec 

(J2000.0) 

^int 

L© 

Distance 

(pc) 

TjoI 

(K) 

DCE 065 


03:28:39.1 

31:06:01.8 

0.02 

250 (50) 

29 (3) 

DCE 004 

L1521F 

04:28:38.9 

26:51:35.6 

0.03 

140 (10) 

20 (3) 

DCE 064 


03:28:32.6 

31:11:05.3 

0.03 

250 (50) 

65 (12) 

DCE 031 

L673-7 

19:21:34.8 

11:21:23.4 

0.04 

300 (100) 

24 (6) 

DCE 001 

IRAM 04191 

04:21:56.9 

15:29:46.0 

0.04 

140 (10) 

27 (3) 

DCE181 


16:28:48.5 

-24:28:38.6 

0.05 

125 (25) 

430 (19) 

DCE 081 


03:30:32.7 

30:26:26.5 

0.06 

250 (50) 

33 (4) 

DCE 025 

L328 

18:16:59.5 

-18:02:30.5 

0.07 

270 (50) 

62 ( 9 ) 

DCE 038 

L1014 

21:24:07.6 

49:59:08.9 

0.09 

250 (50) 

66 (21) 

DCE 185 

IRAS 16253-2429 

16:28:21.6 

-24:34:23.4 

0.09 

125 (25) 

30 (2) 

DCE 109 


03:44:21.4 

31:59:32.6 

0.11 

250 (50) 

345 (4) 

DCE 092 


03:33:14.4 

31:07:10.9 

0.14 

250 (50) 

47 (5) 

DCE 182 


16:27:05.2 

-24:36:29.5 

0.15 

125 (25) 

105 (2) 

DCE 107 


03:44:02.4 

32:02:04.9 

0.15 

250 (50) 

76 (4) 

DCE 063 


03:27:38.3 

30:13:58.8 

0.20 

250 (50) 

199 (3) 

DCE 090 


03:32:29.2 

31:02:40.9 

0.20 

250 (50) 

114 (17) 


Note. — The source informations, Col. (3)-(7), are from [Cuiihari^_et_alJ J2UU2). The order of sources is based on 
their internal luminosities (Col. (5)) and is used in this paper. 


off positions) for each source was 5-15 min. The rms 
noise levels are in a range of 0.06 to 0.13 K at a spec¬ 
tral resolution of 0.078 km s“^. The half power beam 
width [HPBW) is ^67.5" and the pointing accuracy is 
5" rms. The spectra were calibrated in main-beam tem¬ 
perature scale assuming a main-beam efficiency of 0.59 
± 0.07 based on the Saturn measurement (see ARO web 
pagcQ). We used the 3mm (ALMA Type Band 3) receiver 
and the backend Millimeter Autocorrelator with 6.2 kHz 
spectral resolution (^0.019 km s“^), but the spectra were 
smoothed to 0.078 km s“^ for the analysis. 

2.3. SMT observations 

The N 2 H+ (3-2), N 2 D+ (3-2), HCO+ (3-2) and HCN+ 
(3-2) observations were carried out between November 
2012 and February 2014 using the Submillimeter Tele¬ 
scope (SMT) at the ARO, except for N 2 D+ (3-2) obser¬ 
vation for DCE 065 which was carried out in December 
2009. The pointed positions were centered at the infrared 
sources (Table 1) and the observations were performed 
with the position-switching mode. The integration time 
was 6-24 min for N 2 H+ (3-2), 30-42 min for N 2 D+ (3- 
2), 6 min for HCO+ (3-2) and 5 min for HCN (3-2) for 
each source. The rms noise levels are 0.02-0.09 K for 
N 2 H+ (3-2), 0.01-0.04 K for N 2 D+ (3-2), 0.04-0.13 K 
for HCO+ (3-2), and 0.05-0.08 K for HCN (3-2) at a 
spectral resolution of 250 kHz. The HPBWs are 27" for 
N 2 H+ (3-2), 32.6" for N 2 D+ (3-2) and -26" for both 
HCO+ (3-2) and HCN (3-2), and the pointing accuracy 
is 1" rms. We used the ALMA Type 1.3mm receiver 
and the filter-bank backend. The main beam efficiency 
is 0.74 for all linefl The spectral resolution is 250 kHz 
corresponding to 0.27 km s“^ for N 2 H+ (3-2), 0.32 km 
s“^ for N 2 D+ (3-2) and —0.28 km s“^ for both HCO+ 
(3-2) and HCN (3-2). 

3. RESULTS 


^ http://aro.as.arizona.edu/12_obs_manual/appendix_C.htm 
^ http://aro.as.arizona.edu/smt_docs/smt_beam_eff.htm 


All the data were reduced using the CLASS progranfl- 
Eigure 1 shows the spectra of N 2 H+ (1-0), N 2 H+ (3-2), 
N 2 D+ (3-2), HCO+ (3-2) and HCN (3-2) for our targets. 

3.1. N 2 H+ and N 2 D+ 

3.1.1. Spectra 

The N 2 H+ (1-0) lines were detected in all our tar¬ 
gets although some sources have very weak detections. 
We fitted all hyperfine components in N 2 H+ (1-0) by 
considering their relative intensities using our own code. 
The code considers the relative intensities of all compo¬ 
nents similar to the fitting function “METHOD HFS” in 
CLASS. The only difference is that we fit the excitation 
temperature (Tex) directly which gives us a direct esti¬ 
mate of its uncertainty whereas the uncertainty on Tex 
needs to be obtained through error propagation from two 
parameters when using CLASS. Since the “METHOD 
HES” does not provide the covariance of the output pa¬ 
rameters, our fitting process provides a more accurate er¬ 
ror estimate for Tex. Table 2 lists the best fitting results 
including systemic velocities (Vlsr), line widths (AU), 
optical depths (r) and excitation temperatures (Tex), to¬ 
gether with integrated intensities (W). 

The N 2 H''' (3-2) lines were detected toward all the 
targets except DCE 181 and DCE 182 which have weak 
detections in N 2 H''' (1-0). The N 2 D+ (3-2) lines were 
detected toward 11 targets and were undetected in DCE 
181, 025, 038, 182 and 090. The satellite lines are not 
clearly detected with the rms noise a of 0.02-0.09 K for 
N 2 H+ (3-2) and 0.01-0.04 K for N 2 D+ (3-2). This pre¬ 
vents us from deriving the excitation temperatures and 
optical depths (for most of our targets) with HFS fitting. 
However, we can estimate the upper limits of r for each 
source. If we use a threshold 3tT as the upper limit of 
the undetected satellite lines, we obtain a corresponding 
upper limit of r. (The optical depths of N 2 H+ (3-2) of 
DCE 065, 81 and 092 are derived by fitting the hyperfine 
structure with relative intensity.) We thus list the inte¬ 
grated intensities and upper limits of optical depths for 
N 2 D+ and N 2 H+ (3-2) in Table 2. 

^ http://www.iram.fr/IRAMFR/GILDAS/ 








4 


(1-0) _ (3-2) _ (3-2) _ HCO^ (3-2) _ HCN {3-2} 


2 

1 

- DCE 065 

1.0 

0.5 

0.0 


0.6 

0.3 

0.0 


2 

1 

0 


0.6 

0.3 

0.0 

p-'U—^ 'Wlr^nTl/'A^ 

1 

- DCE 004 . 

1.0 

0.5 

0.0 

-^^- 

0.6 

0.3 

0.0 

:-^^-: 

2 

1 

0 

-^^^- 

0.6 

0.3 

0.0 

:-^^^-: 

2 

1 

- DCE 064 ' 

1.0 

0.5 

0.0 


0.6 

0.3 

0.0 

' ' x5 - 



0.6 

0.3 

0.0 


2 

1 

0 

-DCE 0^1 

1.0 

0.5 

0.0 


0.6 

0.3 

0.0 


■ 


0.6 

0.3 

0.0 


2 

1 

- DCE obl ^ ^ 

1.0 

0.5 

0.0 


0.6 

0.3 

0.0 


2 

1 

0 


0.6 

0.3 

0.0 


1 

- DCE I'si j ' x5 - 

1.0 

0.5 

0.0 

-^^- 

xlO - 

% nMVr*’ 

0.6 

0.3 

0.0 

^^^-^ 5 — 

2 

1 

;-^^ 

0.6 

0.3 

0.0 

- -^— r ^ . 

2 

1 

- DCE OSl 

1.0 

0.5 

0.0 


0.6 

0.3 

0.0 


2 

1 

0 

A 

0.6 

0.3 

0.0 


2 2 

h-® 0 

- DCE 0^5 x3 

1.0 

0.5 

0.0 

- ^^—10~ 

0.6 

0.3 

0.0 

^^^-As— 

2 

1 

0 

- ^^^ - 

0.6 

0.3 

0.0 

: - ^^—r ^ 

-AriAp.^J'Al^r1.|All/lr 

2 

DCEOV ' ' x3 

1.0 

0.5 

0.0 


0.6 

0.3 

0.0 


2 

0 

' ' 'x3 

0.6 

0.3 

0.0 

x3 

2 

DCE 1'85 

1.0 

0.5 

0.0 


0.6 

0.3 

0.0 


2 

1 

0 


0.6 

0.3 

0.0 


2 

1 

- DCEibg ' ' 

1.0 

0.5 

0.0 


0.6 

0.3 

0.0 


2 

1 

0 


0.6 

0.3 

0.0 

^TnnjNrir-ij^ 

1 

-DCE 092 ^ 

1.0 

0.5 

0.0 

-^^- 

0.6 

0.3 

0.0 

: -^^-: 

2 

1 

0 

-^^-h- 

0.6 

0.3 

0.0 


2 

1 

0 

DCE 182 ' ' x7 - 

1.0 


0.6 

0.3 

0.0 

: -^^ — Z —7 

x5 

2 

1 

-^^-H- 

0.6 

0.3 

0.0 

-^^- As —^ 

tiApAiiJi f^r^l]n/iA*\- 

2 

DCEIO7 

1.0 

0.5 

0.0 


0.6 

0.3 

0.0 


2 

1 

0 

' A 

0.6 

0.3 

0.0 

' xO.5 ■ 

2 

1 

- DCE 063 

1.0 

0.5 

0.0 


0.6 

0.3 

0.0 

^^^—E—^ 

x5 



0.6 

0.3 

0.0 


1 

- DCE obo ^^5 

1.0 

0.5 

0.0 


0.6 

0.3 

0.0 

^^- As— 

2 

1 


0.6 

0.3 

0.0 



10 -5 0 5 

Velocity-VLSR s'^) 


10 -5 0 5 

Velocity-yLSR (km s'^) 


10 -5 0 5 

Velocity-VLgR (km s'M 


10 -5 0 5 

Velocity-ULSR ^' ^} 


10 -5 0 5 

Velocity-VLSR ^'^} 


Fig. 1.— The spectra of N 2 H’*' (1—0), N 2 H'’" (3—2), N 2 D'^ (3—2), HCO~*" (3-2) and HCN (3-2) for all our LLO targets. The velocities 
are shifted by the Vlsr which is from the N 2 H"^ (1~0) hyperfine fitting (Table 2). The spectra of sources with weak or no detections are 
multiplied by a factor labeled in the upper right corner. 


3.1.2. Two components of DCE 107 

The N 2 H+ (1-0) spectrum of DCE 107 shows a clear 
double peak structure. The SCUBA 850 con- 
tinuum emission im age from the COMPLETE project 
(|Ridge et al. 11200^ has resolved two components with a 
separation of ^ 35" within our observing beam. There¬ 
fore, we suspect that the two velocity components seen in 
the N 2 H+ (1-0) spectrum correspond to the two contin¬ 
uum sources in the 850 /rm image (hereafter DCE 107a 
and DCE 107b). We fitted the spectrum using a su¬ 
perposition of two sets of hyperfine components; each 
set has its own systemic velocity, line width, optical 
depth and excitation temperature (Table 2). The fit¬ 
ting result seems to be reason able by comparing w ith the 
N 2 H+ (1-0) observations from iKirk et al. I (I2007f). Using 
the IRAM 30m telescope {9fwhm 27"). iKirk et al. I 
(I 200 I observed DCE 107 (core 21 in their paper) at a po¬ 
sition ~12" away from our position (the infrared source) 
and ^8" away from DCE107a. They detected only one 
velocity component, supporting our interpretation above, 
and they obtained a line width of 0.65 ± 0.009 km s“^ 
and a systemic velocity of 9.02 ± 0.008 km s“^ which are 
very close to the values we derive for DCE 107a (AU = 
0.61 ± 0.06 km s"! and Ulsr = 8.95 ± 0.02 km s"!). 


3.2. HCO~^ and HCN observations 
3.2.1. spectra 

The HCO+ (3-2) line was detected toward all 16 LLO 
targets and the HCN (3-2) line toward 13 out of the 16 
targets (undetected in DCE 181, 038 and 182). We use 
HCO+ (3-2) and HCN (3-2) as our optically thick trac¬ 
ers to search for infall signatures. The optical depth of 
HCN (3-2) was estimated by fitting the hyperfine com¬ 
ponents. The optical depths of our targets are all larger 
than 5 (except for DCE 025 and DCE 090 which have 
very uncertain values), indicating that the HCN (3-2) 
lines are optically thick toward our targets (Table 3). 
This feature can be easily seen in the spectra since the 
satellite lines (E= 3 —> 3) are significant compared with 
the main components. Note that the satellites E = 2 —>■ 2 
are much weaker or undetected although E = 2 —2 and 
F = 3 —>■ 3 have equal line strengths (anomalies, Daniel 
et al. 2006). 

4. ANALYSIS 

4.1. N 2 H^ Non-LTE RADEX analysis 

We used the non-LTE r adiative transfer code RADEX 
(|van der Tak et al~ll2007f) to model the N 2 H+ (1-0) and 
(3-2) spectra for our targets. The densities of low mass 
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TABLE 2 
LINE PROPERTIES 





N 2 H+ (1-0) 



N 2 H+ (3-2) 

N 2 D+ (3-2) 


-^2-ff j = i_o 

Vlsr 

Ay (FWHM) 

"Ttot^ 

q-i b 
-L ex 

W 4- 

^2^j = 3_2 

'Ttot'^ 

W -4- 

-^2^j=3_2 

rtot*^ 

Source 

[K km s-l] 

[km s-l] 

[km s~^] 


[K] 

[K km s-l] 


[K km s-l] 


DCE 065 

3.58 (0.07) 

7.02 (0.01) 

0.46 (0.02) 

5.00 (1.13) 

4.95 (0.30) 

0.69 (0.04) 

2.84 (1.10) 

0.34 (0.04) 

<13.10 

DCE 004 

3.19 (0.08) 

6.44 ( 0 . 01 ) 

0.29 (0.01) 

14.05 (2.22) 

4.34 (0.09) 

0.52 (0.06) 

<12.30 

0.19 (0.02) 

<6.60 

DCE 064 

2.31 (0.08) 

7.19 (0.01) 

0.35 (0.02) 

2.94 (1.27) 

5.52 (0.87) 

0.56 (0.06) 

<18.70 

0.06 (0.02) 

<19.90 

DCE 031 

1.62 (0.08) 

7.15 (0.01) 

0.51 (0.03) 

2.07 (1.30) 

4.60 (0.91) 

0.55 (0.04) 

<8.20 

0.11 (0.01) 

<14.00 

DCE 001 

2.28 (0.07) 

6.63 (0.02) 

0.62 (0.04) 

3.05 (1.49) 

4.39 (0.56) 

1.29 (0.03) 

<2.40 

0.29 (0.01) 

<4.10 

DCE 181 

0.33 (0.06) 

3.67 (0.01) 

0.26 (0.03) 

0.10 (0.07) 

14.48 (18.46) 

<0.02^ 

- 

<0.02^ 

- 

DCE 081 

2.62 (0.07) 

6.05 (0.01) 

0.40 (0.02) 

5.79 (1.44) 

4.43 (0.24) 

0.56 (0.03) 

8.68 (2.28) 

0.14 (0.01) 

<7.30 

DCE 025 

0.95 (0.07) 

6.61 (0.01) 

0.40 (0.03) 

4.43 (2.12) 

3.51 (0.24) 

0.11 (0.02) 

<19.90 

<0.02=* 

- 

DCE 038 

0.40 (0.05) 

4.24 (0.03) 

0.44 (0.07) 

1.95 (4.19) 

3.40 (1.16) 

0.07 (0.02) 

<19.90 

<0.03'i 

- 

DCE 185 

2.10 (0.05) 

4.04 (0.004) 

0.23 (0.01) 

5.76 (1.51) 

4.93 (0.34) 

0.73 (0.03) 

<19.50 

0.21 (0.01) 

<4.90 

DCE 109 

2.07 (0.07) 

9.12 (0.01) 

0.45 (0.02) 

0.87 (1.31) 

8.49 (7.59) 

0.32 (0.04) 

<12.10 

0.14 (0.02) 

<7.50 

DCE 092 

6.46 (0.08) 

6.27 ( 0 . 01 ) 

0.66 (0.02) 

6.89 (0.88) 

5.12 (0.14) 

1.48 (0.05) 

5.06 (1.33) 

0.59 (0.02) 

<0.90 

DCE 182 

0.32 (0.05) 

4.84 (0.03) 

0.37 (0.07) 

5.55 (6.20) 

2.97 (0.16) 

<0.03=^ 

- 

<0.02^* 

- 

DCE 107a 

1.34 (0.85) 

8.95 (0.02) 

0.62 (0.06) 

1.57 (0.46) 

4.61 (0.76) 

1.04 (0.04) 

<3.90 

0.23 (0.02) 

<4.50 

DCE 107b 

0.87 (1.66) 

8.40 (0.01) 

0.33 (0.03) 

1.08 (0.87) 

6.30 (4.45) 

1.04 (0.04) 

<3.90 

0.23 (0.02) 

<4.50 

DCE 063 

2.59 (0.08) 

4.71 (0.02) 

0.65 (0.04) 

1.35 (1.31) 

5.92 (2.53) 

0.84 (0.04) 

<4.50 

0.13 (0.02) 

<19.90 

DCE 090 

0.37 (0.04) 

6.62 (0.01) 

0.29 (0.03) 

0.10 (1.79) 

15.15 (325.44) 

0.21 (0.03) 

<19.90 

<0.02^ 

- 


Note. — The values in parenthesis represent 1<7 uncertainties. 


T represents the total opacity which is the sum of the optical depth of all hyperfine components. The fitting limitation of t is set to be t > 0.1. 

The excitation temperatures are computed by assuming a beam filling factor of 1. 

The upper limits of the total opacity t are estimated based on the peak intensities and the rms noise levels since most of the satellite lines are undetected. 
The integrated intensity upper limits is estimated using a peak temperature detection limit of Scr, and assuming the line width is the same as N 2 H'^ (1—0). 


TABLE 3 

Line properties of HCO+ (3-2) and HCN (3-2) 


Source 

HCO+ (3-2) 

’^nco+_3_2 
[K km s-i] 



HCN (3-2) 




Vlsr 

[km s“^[ 

AV 

[km s“^[ 

Ttot^ 

rn C 

J- ex 

[K] 

WncJVj=3_2 
[K km s-^[ 

DCE 

065 

5.96 (0.03) 

6.92 (0.01) 

0.27 (0.03) 

28.04 (9.85) 

3.98 (0.10) 

2.37 

(0.04) 

DCE 

004 

6.52 (0.07) 

6.30 (0.02) 

0.37 (0.04) 

43.28 (18.33) 

3.75 (0.09) 

1.99 

(0.07) 

DCE 

064 

7.25 (0.05) 

7.05 (0.05) 

0.47 (0.09) 

13.19 (7.46) 

3.50 (0.09) 

1.43 

(0.06) 

DCE 

031 

4.37 (0.05) 

7.17 (0.05) 

0.72 (0.13) 

5.16 (3.04) 

3.82 (0.11) 

1.81 

(0.07) 

DCE 

001 

8.23 (0.07) 

6.66 (0.05) 

0.56 (0.07) 

19.64 (8.84) 

3.54 (0.08) 

1.50 

(0.07) 

DCE 

181 

0.19 (0.02) 

- 

- 

- 

- 


- 

DCE 

081 

4.00 (0.04) 

6.01 (0.05) 

0.57 (0.07) 

29.90 (13.24) 

3.54 (0.09) 

1.83 

(0.05) 

DCE 

025 

2.93 (0.04) 

6.71 (0.21) 

1.30 b-05) 

1.21 (7.07) 

3.34 (1-68) 

0.60 

(0.06) 

DCE 

038 

1.36 (0.03) 

- 

- 

- 

- 


- 

DCE 

185 

6.10 (0.07) 

4.05 (0.03) 

0.32 (0.05) 

36.18 (15.39) 

3.87 (0.10) 

2.01 

(0.09) 

DCE 

109 

5.63 (0.03) 

9.13 (0.03) 

0.54 (0.05) 

17.57 (5.14) 

3.97 (0.08) 

2.73 

(0.04) 

DCE 

092 

10.06 (0.10) 

5.82 (0.04) 

0.44 (0.07) 

12.74 (6.40) 

3.76 (0.10) 

2.36 

(0.11) 

DCE 

182 

7.63 (0.06) 

- 

- 

- 

- 


- 

DCE 

107 

19.73 (0.07) 

8.91 (0.03) 

0.73 (0.06) 

8.20 (2.29) 

4.37 (0.07) 

4.38 

(0.10) 

DCE 

063 

8.18 (0.05) 

4.45 (0.05) 

0.67 (0.08) 

14.09 (5.77) 

3.51 (0.07) 

1.73 

(0.05) 

DCE 

090 

3.03 (0.04) 

6.42 (0.19) 

0.29 (0.55) 

0.10 (3.50) 

9.32 (174.78) 

0.67 

(0.04) 


Note. — The values in parenthesis represent Itr uncertainties. 

^ T represents the total opacity which is the sum of the optical depth of all hyperfine components. The fitting limitation of 
T is set to be t > 0.1. 

The excitation temperatures are computed by assuming a beam filling factor of 1. 


protostellar cores can be lower than the critical densi¬ 
ties of N 2 H+ (1-0) and (3-2). This means that the local 
thermal equilibrium (LTE) assumption is not appropri¬ 
ate for deriving the true physical properties. RADEX 
provides fast non-LTE analysis for molecular line spec¬ 
tra in a uniform medium, which involves collisional and 
radiative processes. The N 2 H+ line data (energy lev¬ 
els, statistical weights, Einstein A -coefficients a n d col ¬ 
lisio nal rate coeffic i ents) are from iDaniel et al. I l)2005f) 
and lSchoier et alTl (|2005D . and were taken from the Lei¬ 


den Atomic and Molecular Database (LAMDA). 

We used RADEX to derive the cores’ kinetic tempera¬ 
tures (Tkin), N 2 H+ column densities (N(N 2 H+)), H 2 den¬ 
sities (riHa) s-iid source sizes. The former three param¬ 
eters are related to the cores’ evolutionary state and/or 
may decide the mass of stars formed in the future. The 
core sizes, together with the beam sizes, yield the beam 
filling factors that determine the scale of the spectra. 
Since the beam sizes of our N 2 H+ (1-0) and (3-2) obser¬ 
vations are different, the beam filling factors affect the 
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intensity ratios of N 2 H+ (l-0)/(3-2) as well as our de¬ 
rived physical parameters. 

We constructed a model grid for each source to de¬ 
rive the physical parameters. The grid is composed with 
three variances, kinetic temperature, N 2 H“*' column den¬ 
sity and H 2 density, while the line width is fixed to that 
from the hyperfine fitting. The step sizes of the grid are 
0.5 K for Tkin, 0.025 and 0.05 in decimal log scale for 
N(N 2 H''') and riHa, respectively. Each cell in the grid 
contains the synthetic spectra of N 2 H+ (1-0) and (3- 
2). The synthetic spectra are built-up based on RADEX 
results which include the excitation temperatures (Tex) 
and optical depths (r) for all hyperfine components. The 
spectra are constructed using the equation 

Tmb{v) = vI,(^l®i^-J(Tbg))(l-exp(-^r,(u))), 

( 1 ) 

where 'f' is the beam hlling factor, Tbg is the cos¬ 
mic background temperature (2.73 K), and J{T) = 
[hi//k]/[exp{hi//kT) — 1]. To approach a reasonable in¬ 
tensity at the hyperfine component overlapping region, 

we used a weighting emissivity term, ^ . This 

conception is similar to t he av eraged excitation tem¬ 
perature in iDaniel et al. I (|2006ll for discussing the line 
anomalies {i.e. Tex are different between the hyper¬ 
fine components.). The weighting is based on the op¬ 
tical depths contributed from all components at a spe¬ 
cific velocity. X]a(^^) is the superposition of optical 
depths contributed from all components at such veloc¬ 
ity. Note that RADEX does not consider the line over¬ 
lap, whi ch may make our approach uncertain at suc h 
regions (jDaniel et al. 1 120061 : | van der Tak et al. ]| 2003 ). 
The beam filling factor {i.e. source size) is the free pa¬ 
rameter determining the scale of the synthetic spectra. 
As a result, we obtained a distribution (see Appendix 
and Figure Al) and the corresponding best fitting core 
size distribution in the grid for each source. Figure A2 
shows an example (DCE 081) of the best fitting results 
using our RADEX synthetic spectra compared with us¬ 
ing HFS fitting. The RADEX synthetic spectra seem to 
fit the hyperfine anomalies better than HFS fitting for 
the 121-011 and 122-011 components. For the 112-012 
component, both htting processes can not ht the anoma¬ 
lies well. For most of our sources, no significant difference 
is found between the two fitting processes. The best ht¬ 
ting results are listed in Table 4 for the 13 targets with 
both N 2 H+ (1-0) and (3-2) detections and only one ve¬ 
locity component. T he errors are dehii ed with a 68.3% 
conhdence level area (|Press et al. 1119921) for four free pa¬ 
rameters (Tkin, N(N 2 H+), riHa and core sizes). 

Out of the 13 targets, Tkin and are constrained in 
5 sources but not in 8 sources (which have a conhdence 
interval exceeding the boundary of our model grid, see 
Figure Al). For the 5 sources with constraints, the av¬ 
erage error is '^3.1 K for Tkin and a multiplicative factor 
of ~0.36 for riHs (steps by log scale). DCE 064 has the 
most uncertain Tkin which is IS.blg'g K and DCE 065 

has the most uncertain uhs which is 10 -o '! cm”'^. For 
the 8 unconstrained sources, we list the parameters from 
the best-ht model and upper or lower limits of these pa¬ 
rameters (Table 4). 


TABLE 4 

RADEX FITTING RESULTS 



Tkin 


N(N2H+) 

FWHM 

Source 

[K] 

logic [cm-3] 

logic [cm-2] 

[arcsec] 

DCE065 

7 c:+2.5 

^•-*--0.4 

13.10tn2 

58tll 

DCE004 

gO+X5 

K Q + 00 

O.o_o.4 

1 Q oc+0.15 
iO.oo_0.30 

45lf 

DCE064 



13.021®+! 

42 I!® 

DCE031 

OK C + OO 

zo.o_23 Q 

5 . 0 I?,® 

13.23l®i! 

27ll 

DCEOOl 


7 q+0.2 

' •'^-0.7 

1 0 -| 7+0.06 
10.1 f _o.07 

4211 

DCE081 


c 9+0.3 
O.Z_o 2 

1 09 + 0.13 

10.0Z_Q Qg 

3411 

DCE025 


c ^+00 

1 9 Q9+0.31 

iz.oz_Q g2 

4ol!! 

DCE038 

9.0l^ 

5 . 0 +“ 

12.951®;®® 

24l“ 

DCE185 

6.5+?® 

■7 q+00 

' •'^-1.3 

12.781®;!! 

77+18 

' ' -34 

DCE109 

75 +K .5 

6-otro 

12.631®;!® 

85lS 

DCE092 

7 c+O-S 

' -^-0.5 

6.3l°;? 

1 0 07+0.06 
lO.Of _Q Q2 

65+! 

DCE063 

33.0t^,5 


1 90 + 0.09 

io.zo_0.33 

35li® 

DCE090 

30.5l^,o 

52+0^8 

13.281®;!® 

111 ® 


Note. — The parameters of the best-fit model are listed in the table. The 
errors presented here are with 68.3% confidence level. The 00 error indicates 
that the parameter is not constrained at the upper or lower end. 


We here compare the H 2 densities of our fitting results 
with the literature for several well-studied objects. Our 
derived densities are quite consiste nt with that o f DCE 
004 fm 521F) and 038 (L1014) in ICra.nsi et al.~l ([ 200 I 
l2005ajl . Considering the 1.2 mm continuum maps and a 
density profile, = no/(l-l-r/ro)“, they found no = 
10® cm-3 and a - 2 for DCE 004, and no = 2.5 x 10® 
cm“® and a ~ 2.7 for DCE 038 (both with a “flat” re¬ 
gion, ro ^ 20"). Our derived densities (6.3 x 10® cm“® 
within r < 23" for DCE 004, and 10® cm“® within r < 
12" for DCE 038) are very close to the results from the lit¬ 
erature (but DCE 038 is not well constrained in our sam¬ 
ple) . However, the particularly high densities of DCE 001 
(IRAM 04191) and DCE 185 ( IRAS 16253) are very li kely 
overestimated. For DCE 001. iBelloche et al~l (l2002t) de¬ 
rived a density of ^ 2 x 10® cm“® at a radius of 20", based 
on the dust continuum emission and the radiative trans¬ 
fer modeling of CS transitions. Our estimated density, 
2 X 10^ cm“® within a radius of ^21" is obviously too 
high. DCE 001 has the highest N 2 H+ (3-2)/(l-0) inten¬ 
sity ratio among our targets (Table 2). In addition, our 
observed N 2 H+ (1-0) peak int ensity is lower than t hat 
from the IRAM 30m telescope (jBelloche et al. l[2002h by 
a factor of ^3. The optical dept h derived from the HFS 
fitting is also lower than that in IBelloche et al.~l l)2002[) 
by a factor of ^2. We therefore speculate that our obser¬ 
vation with a large beam {HPBW = 67.5") is seriously 
contaminated by diffuse regions whereas the IRAM 30m 
observation {HPBW ^27") is not. Our N 2 H+ (3-2) ob¬ 
servations have the same angular resolution as the IRAM 
30m observations of N 2 H+ (1-0), so they likely trace the 
same region. The contamination of our N 2 H+ (1-0) data 
by more diffuse material likely leads to an overestimate 
of the N 2 H+ (3-2)/(l“0) intensity ratio that is used in 
our analysis to derive constraints on the density. Tak¬ 
ing into account beam filling factor effects as we did in 
our analysis cannot solve this problem. The case of DCE 
185 is very similar. The high N 2 H+ (3-2)/(l-0) inten¬ 
sity ratio is likely due to the contamination of our 1-0 
spectrum by more diffuse regions; the high density (2 x 
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10^ cm“^) within a radius of 39" are unlikely to be true 
since the envelope mass of DCE 185 is small (0.15-1.0 
Mq Barsony et al. 2010, Stanke et al. 2006, Enoch et 
al. 2008 and Tobin et al. 2012). Some sources with high 
N 2 H+ (3-2)/(l-0) intensity ratios may also suffer from 
this problem such as DCE 031 and 063 (intensity ratios 
3). As a result, these two sources have the highest 
kinetic temperatures among our targets for explaining 
their high intensity ratios in our RADEX fittings. The 
kinetic temperatures could thus be overestimated. It is 
likely that sources with relatively low N 2 H“'' (3-2)/(l-0) 
intensity ratios (and/or high opacities) are less affected 
by the contamination of diffuse regions. N 2 H+ (3-2) and 
(1-0) observations with a similar beam size are required 
for deriving the densities of DCE 001, 185, 031, and 063 
in a more reliable way. 


TABLE 5 

N 2 D+/N 2 H+ COLUMN DENSITY RATIOS 


Source 

N(N2H+) 
10^^ cm“^ 

N(N2H+) 
I 0 I 2 cm-2 

N(N2D+) 

lO^l cm-2 

N(N2D+)/N(N2H+) 

DCE 

065 

16.73 

(0.96) 

27.54 

(1.39) 

23.91 

(2.81) 

0.087 

(0.011) 

DCE 

004 

33.94 

(2.24) 

88.26 

(10.88) 

23.96 

(2.92) 

0.027 

(0.005) 

DCE 

064 

14.15 

(1.69) 

27.95 

(3.20) 

5.03 

(1.46) 

0.018 

(0.006) 

DCE 

031 

21.38 

(3.38) 

28.53 

(2.18) 

11.43 

(1.63) 

0.040 

(0.006) 

DCE 

001 

8.53 

(1.41) 

24.96 

( 0 . 54 ) 

13.97 

(0.40) 

0.056 

(0.002) 

DCE 

081 

28.96 

(2.43) 

54.71 

(0.27) 

15.50 

(0.82) 

0.028 

(0.002) 

DCE 

025 

6.02 

(1.85) 

14.04 

(2.76) 


- 

- 


DCE 

185 

6.72 

(0.48) 

11.83 

(0.48) 

7.55 

(0.46) 

0.064 

(0.005) 

DCE 

109 

5.04 

(0.60) 

6.23 

(0.69) 

10.58 

(1.21) 

0.170 

(0.027) 

DCE 

092 

25.93 

(1.06) 

50.59 

(1.68) 

30.68 

(0.86) 

0.061 

(0.003) 

DCE 

063 

16.75 

(2.30) 

24.12 

(1.10) 

9.41 

(1.34) 

0.039 

(0.006) 


Note. — DCE 038 and DCE 090 column densities are not derived because their N 2 H^” 101-012 
component is not detected (see text for details). 


4.2. N2D'^/N2H^ column density ratio 

The N 2 D+ column density is difficult to derive with 
only the J = 3 —2 transition because we cannot derive 
the excitation temperature (Tex) and optical depth (r) 
which are crucial for determining the column density (see 
§3.1.1). In addition, the excitation temperatures could 
be very different between different transitions under non- 
LTE conditions (see Fig. 17 in Daniel et al. 2007). Us¬ 
ing N 2 H+ as an example, the excitation temperatures of 
N 2 H+ (3-2) are lower than that of N 2 H+ (1-0) by 0.3 - 
4 K in our RADEX best-fit models (see §4.1). 

With only one transition, we have to perform an LTE 
analysis to derive the column density of N 2 D+. How¬ 
ever, a reasonable N 2 D+/N 2 H+ column density ratio 
(N( N 2 D~*')/N(N 2 H + )) ca n still be obtained considering 
that iDaniel et al. I (|2007f) found very similar excitation 
temperatures for N 2 D+ and N 2 H''' in a given transition 
(J = l—>-0orJ = 2—>-lorJ = 3—>-2, see their 
Fig. 17). For the J = 3 ^ 2 transition, the largest dif¬ 
ference in excitation temperature is < 0.5 K at a H 2 
density of ^ 10^ cm“^. We therefore assume that the 
N(N 2 D+)/N(N 2 H+) ratio derived using the same transi¬ 
tion (J = 3 —>■ 2) in both molecules will be very close to 
the actual ratio, even when deriving each column den¬ 
sity with LTE equations. We use the excitation temper¬ 
atures of N 2 H+ (3-2) obtained from the RADEX best-fit 
model for both N 2 H+ and N 2 D+ (3-2). (Since each hy- 
perfine component of N 2 H+ (3-2) has a different Tex, 
we use the averaged excitation temperature as defined 

in IDaniel et al. I (| 2006 ll . Tave = 

The opacity of N 2 D+ (3-2) is then calculated by fitting 
the hyperfine components with Tex fixed to the value 
derived for N 2 H+ (3-2). Note the N 2 D+ (3-2) spectra 
were scaled with the beam filing factors derived from 
source sizes. The column de n sities are derived using Eq. 
1761 in iM angnm fc Shirlevi (|201,3H and equation A4 in 
iCaselli et al. ()2002| 1. 


Afthick 


= iV, 


thin 


r 

1 — exp{—T) 


( 2 ) 


and 

^thin ^ SttIU gi _ 1 

A3A [J(Tex) - J(Tbg)] 

1 Qrot(Tex) 

1 - exp(-/nz//cTex) gi exp(-Ti/A:Tex) 

(3) 

where r is the total opacity of the multiplet, W is the 
total integrated intensity of the multiplet (scaled by the 
beam filling factor -L see §4.1), A is 

the Einstein coefficient, gi and are the statistical 
weights of the lower and upper levels, Qrot is the rota¬ 
tional partition function under LTE conditions, and Ei 
is the energy of the lower level. The column densities 
are listed in Col. (3) and Col. (4) of Table 5 and the 
N(N 2 D+)/N(N 2 H+) ratios are listed in Col. (5). The 
uncertainties reported in Col. 5 take into account the 
uncertainties on the measured integrated intensities but 
neglect the uncertainty on Tex because it does not signif¬ 
icantly affect the N(N 2 D+)/N(N 2 H+) ratio. For a typ¬ 
ical excitation temperature of ~ 5.4 K in our sample, 
the N(N 2 D+)/N(N 2 H+) ratio will only increase by ^ 7% 
and decrease by ~ 9% if the excitation temperature in¬ 
creases by 1 K and decreases by 1 K, respectively. The 
effect is even smaller at a higher excitation temperature; 
the N(N 2 D''')/N(N 2 H+) ratio will increase by ^20% if 
the excitation temperature increases by 5 K. Another 
noteworthy uncertainty could come from the unresolved 
hyperfine structure of J = 3^-2. We note also that 
the opacity correction applied in Eq. (2) slightly over¬ 
estimates the column density because the total opacity 
of the multiplet, t, is higher than the opacities of the 
individual components. 

We also derive the N 2 H+ column densities from the J 
= 1 —0 transition with Eqs. (2) and (3), using only the 
101-012 component that is not affected by other hyper¬ 
fine components (see Col. (2) of Table 5). The N 2 H+ 
column densities derived from N 2 H+ (1-0) are slightly 
larger than that from RADEX (the best-fit models) by 
factors of 1.0 to 1.5 except for DCE 001 (0.6) and DCE 
025 (0.9). The average of the factors is 1.2 with a stan¬ 
dard deviation of 0.25. This result makes sense because 
LTE analysis will overestimate the N 2 H+ populations in 
J > 1 levels when the collision rate is insufficient to pump 

























N 2 H'^ with low volume densities (non-LTE). The column 
densities derived from N 2 H+ (3-2) (Col. (3) of Table 5) 
are also larger than that from the RADEX best-fit model 
(Col. (4) of Table 4) by factors of 1.4 - 3.9, with an av¬ 
erage value of 2.1 and a standard deviation of 0.6. The 
column densities of N 2 D+ listed in Col. (4) of Table 5 are 
thus likely overestimated by similar factors. However, as 
explained above, the N 2 D+/N 2 H+ column density ratios 
in Col. (5) should be reliable. 


4.3. Asymmetry parameters 

To probe the infall motion of our tar get cores, we used 
the as ymmetry parameter defined by Mardones et al.1 

c- ^thick "^thin / ., 

5v = ---, (4) 

Authin 


where Uthick and Uthin are the velocities of the optically 
thick line and optically thin line, respectively, and Authin 
is the line width of the optically thin line. We used 
N 2 H+ ( 1 - 0 ) as our optically thin tracer, and Uthin and 
Authin are derived from HES fitting (Table 2). The HFS 
fitting provides the line width that corresponds to the 
width of the opacity profile. It should represent the ve¬ 
locity dispersion well and can thus be used as a proxy 
for Authin- Eleven out of our sixteen sources have op¬ 
tical depths lower than 1.3 for the strongest hyperfine 
component (123-012, rtot x 7/27), suggesting that the 
systemic velocity derived from the HES fitting should not 
be significantly affected by optical depth effects and can 
be used as a proxy for Vthin- For the remaining sources, 
three (DCE 081, 185, and 092) have optical depths lower 
than 1.3 for the isolated hyperfine component (101-012, 
Ttot X 3/27). The systemic velocity derived from the 
HFS fitting matches well the centroid velocity of the iso¬ 
lated component for these three sources and can thus 
also be used as a proxy for Uthin- For DCE 182, although 
the N 2 H+ ( 1 - 0 ) emission is very weak and the opacity 
is not constrained, its emission is most likely optically 
thin. The remaining source, DCE 004, is optically thin 
(r ^ 0.5) in the faintest hyperfine component (110-011, 
Ttot X 1/27). In this source, this component is well 
separated from the other hyperfine components and a 
simple Gaussian fitting yields a systemic velocity consis¬ 
tent with the one derived from the HFS fitting. HCO+ 
(3-2) and HCN (3-2) are used as optically thick tracers. 
Dhick corresponds to the peak velocity and we assume 
a 2a uncertainty of 0.14 km s“^, which is equal to half 
of the channel width. The asymmetry parameters cal¬ 
culated from HCO+ (3-2) and HCN (3-2) are listed in 
Table 6 . 

We use |(5u| > 30 - 5 .;, (S/N > 3 ) to define the significant 
detections of blueshifts (infall motion) or redshifts (out¬ 
ward motion). These sources are separated into three 
groups, which are 1) Sv < —3asv- clear indication of in¬ 
fall, 2) Sv > 3aSv'- redshifted optically thick line, and 
3) 3aSv > Sv > —3asv- no clear indication of systematic 
motion. With the numbers of sources in the three groups, 
we cal culated the blue excess defined bv iMardones et alTI 

eng), 

ui -^blue -^red /r\ 

biue excess =-—-, (5) 


TABLE 6 
Infall signatures 


Source 


5v 


Hn 

[km s“^] 

HCO+ 

HCN 

consistency 

DCE065 

-0.06 (0.15) 

-0.34 (0.16) 

YES 

0.73 (0.17) 

DCE004 

-0.78 (0.24) 

-1.23 (0.25) 

YES 

0.14 (0.02) 

DCE064 

-0.80 (0.20) 

-0.37 (0.20) 

YES 

-0.62 (0.15) 

DCE031 

0.53 (0.14) 

0.63 (0.14) 

YES 

- 

DCEOOl 

-0.36 (0.12) 

-0.57 (0.12) 

YES 

0.07 (0.03) 

DCE181 

0.43 (0.28) 

- 

- 

- 

DCE081 

0.82 (0.18) 

-0.21 (0.18) 

NO 

-0.13 (0.03) 

DCE025 

1.32 (0.20) 

1.43 (0.20) 

YES 

-0.08 (0.04) 

DCE038 

-0.39 (0.18) 

- 

- 

0.41 (0.18) 

DCE185 

1.31 (0.31) 

-1.00 (0.31) 

NO 

- 

DCE109 

0.16 (0.16) 

-0.13 (0.16) 

NO 

-1.04 (0.15) 

DCE092 

-0.39 (0.11) 

-0.59 (0.11) 

YES 

0.58 (0.06) 

DCE182 

-0.58 (0.23) 

- 

- 

0.59 (0.29) 

DCE107 

0.36 (0.13) 

0.15 (0.12) 

YES 

-0.46 (0.08) 

DCE063 

-0.03 (0.11) 

0.20 (0.11) 

NO 

0.72 (0.22) 

DCE090 

-0.08 (0.25) 

-0.54 (0.26) 

YES 

0.26 (0.05) 


Note. — The values in parenthesis represent la uncertainties. 


The blue excess is 0.07 (excluding DCE 185, see below) 
and 0.15 for HCO+ (3-2) and HCN (3-2), respectively. 
These results are not very different f rom that the re¬ 
sults o f previous works, i.e. 0.25 f r om iMardones et al. I 
(I 1997 II. 0.28 from lGregers en et al. I (|200r)H and 0.05 from 
lEmnrechtinger et al. I (|2009H . We also use the binomial 
test to deter mine if the numbe r of blueshifted sources 
is significant (|Rvgl et al. 11201^ . The survival functions 
are 0.50 for HCO+ (3-2) and 0.34 for HCN (3-2), sug¬ 
gesting that the number of blueshifted detections is not 
significant in our sample for both HCO’*' (3-2) and HCN 
(3-2). 


4.4. Two-layer model fitting 

We fi tted the HCO+ (3~ 2) spectra with the two-layer 
model (|Mvers et al. 1 119961 1 to derive the infall veloci- 
tie s of our tygets. We used a simple version defined 
bv iLee et al. I pOOlI l 

Tmb = J(Td(l - e-"’')e-"^ + J(rbg)(e-"^-"^ - 

( 6 ) 

where Tbg is the cosmic background temperature (2.73 
K), Tv is the excitation temperature of the rear layer, 
and Tf and Tv are the optical depths of the front and rear 
layers, respectively. The optical depths are expressed 
as Tf = roexp[—(u — V/)^/2 (t^] and Tv = roexp[—(u -|- 
K)^/2cr^], where tq is the peak optical depth, Vf and W 
are the systemic velocities of the front and rear layers, re¬ 
spectively, and a is the velocity dispersion. Considering 
the relative velocity of the two layers as the infall velocity 
{i.e. Vf = Vr = ldn/2), we fitted the spectra with eq. (6) 
to obtain the infall velocities (Figure 2). A positive infall 
velocity represents an inward motion in the core while a 
negative one indicates an outward motion. The spectral 
resolution (0.28 km s“^) is sometimes insufficient to re¬ 
solve the double-peak features and/or asymmetry struc¬ 
tures. We therefore only list the fitting results with infall 
or outward velocities above the uncertainties (Table 6). 
As a result, the infall velocities are positive in 8 sources, 
negative in 5 sources and unconstrained in 3 sources. 


to study the statistical properties of our LLO sample. 


5. DISCUSSION 
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Fig. 2.— Two-layer radiative transfer model fits to the HCO'*" (3—2) spectra for the 16 targets. The fitting results are shown in blue line 
for Vin > 0, red line for Vin < 0 and green for unconstrained Vin- The vertical dashed line indicates the Flsr from the two-layer model 
fitting and the vertical dotted line indicates the Vlspt from the N 2 H+ hyperfine fitting. 


5.1. The physical conditions from Non-LTE RADEX 

analysis 

Table 4 shows the physical parameters derived with 
RADEX. The kinetic temperature (Tkin), N2H+ column 
density (N(N2H+)), H2 density (riHa) and source sizes of 
five cores are constrained in the RADEX analysis. Eight 
sources have no constraints in their kinetic temperatures 
and H2 densities; the data are unable to discriminate 
models between a wide variety. Additional spectral lines 
or higher angular resolution data would be required to 
discriminate the models for these sources. Data with 
higher sensitivity could also break the degeneracy (see 
Appendix). The variances of Tkin and nua from the fit¬ 
ting models are always anti-correlated; the fitting-models 
can either have high Tkin with low riHa "^kin with 

high riHa- With the currently available data, we dis¬ 
cuss the physical properties with the best-fit models from 
both the 5 constrained sources and the 8 unconstrained 
sources. Our LLO targets have Tkm in a range between 
5.0 K and 33 K, and a median of 8.0 K. (The correspond¬ 
ing thermal line widths are 0.09, 0.24 and 0.12 km s“^ for 
N2H+.) The H2 densities are populated in a range be¬ 
tween 10® cm“® and 2x10® cm“® with a median of 4x10® 
cm“® (excluding DCE 001 and 185 which have uniquely 
high riHa = 2x10^ cm“® in our sample, see §4.1). Com¬ 
paring with the Class 0 sample ( Tkin = 9 - 15 K and 
rtHo = 5x10® - 10^ cm“®) from lEmorechtinger et al~~l 
(|2009H . LLOs have relatively lower riHa and a wider Tkin 
distribution; the standard deviation of Tkin for our LLO 


sample is 9.3 which is much larger than that of the 
Class 0 sample, ^ 1.7. The derived H2 densities of LLOs 
also spread in a larger range (~1.3 magnitude) than the 
Class 0 objects even after excluding DCE 001 and 185. 
We think this result comes from a wide variety of natures 
of our LLO targets which include some Class I objects 
and are located in different molec ular clouds or isolated 
cores while the Class 0 sources in lEmorechtinger et al. I 
(I 2 OO 9 I) are all in the P e rseus molecular cloud. Note that 
lEmorechtinger et al. I (I2009D derived riHa using N2D+ 
(2-1) and (3-2) lines, and the kinetic temperatures were 
derived bv iHatchell I (|2003D and iJiiina et al. I (jl999D with 
NH3 lines. 

A caveat in the H2 density comparison is that the 
beam size of our dat a is m uch larger than that from 
lEmorechtinger et al. I 1)20091) . The larger beam could 
make the derived H2 density lower if it is diluted by a 
large area. To fairly compare the H2 densities, we scale 
our derived nuo to the observing beam size o f N2D+ (2- 
1) (16.3") from lEmorechtinger et al71 (I2009D . Assuming 
a density profile, p(r) oc r“^, we calculate an average den¬ 
sity, Pave(<i') oc r“^, within a radius r. We can therefore 
obtain an expected H2 density within 16.3" area based 
on our derived H2 density and source size. The scaled 
H2 density is populated in a range between 7x10^ cm“® 
and 3x10^ cm“® with a median of 2x10® cm“®. Our 
derived H2 densities therefore have a wider distribution 
compared with the Class 0 sample. Excluding DCE 001 
and 185 (see §4.1), 3 LLOs have H2 density higher than 
Class 0 sources and 8 have lower H2 density. There- 
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Fig. 3.— Comparison of N 2 D+/N 2 H+ column density ratio with Ti;iu and The Class 0 data are from lEmprechtinger et al. 1 1200911 
(green) and IFriesen et al. I 120131 (yellow). The error bar indicates Icr uncertainty in our data (red). The green area represents the range 
of H 2 densities of the Class 0 objects from lEmprechtinger et al. 1 1200911 . 
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Fig. 4.— Populations of line width (a) and N 2 D'*'/N 2 H'*" column density ratio (b) for Class 0 p r otostars (green) , Starl ess co res (blue) and 
LLOs (red). The Class 0 protostar and starless core data are from lEmprechtinger et al. I 12009l ~l. [Friesen et al. I 12013 ~) and ICrapsi et al~l 
12005bll . 


fore, LLOs are likely to have relatively low H2 densities 
compared with Class 0 sources after correcting for the 
different beam sizes. 

5.2. N(N 2 D'^)/N(N 2 H^) ratio as an age indicator 

To investigate the relation between chemistry and 
physical conditions in these cores, we compare the 
N(N2D+)/N(N2H+) ratio with the derived kinetic tem¬ 
perature (or gas temperature) and H2 density (Fig- 
ure 3). Figure 3a includes the Class 0 objec t s from 
lEmprechtinger et al. I (j2009[) and iFriesen et al. I (j201,‘lD . 
The kinetic temperatures are derived using N2H''" J = 3- 
2/J = 1-0 in our study, and usin g NH.^j (1, 1)/ (2 2) i n 
lEmprechtinger et al. I (|2009D and IFriesen et al. 1 (|2013D . 
The sources with the highest N(N2D+)/N(N2H+) ratios 


(> 0.05) all have low kinetic temperatures (< 10 K). We 
use Pearson’s r c orrelation test and Sp earmans’s p rank 
correlation test (jConover et al. 1 119991) to evaluate the 
significances of correlations between N(N2D+)/N(N2H''') 
ratios and kinetic temperatures for our LLO sample. The 
Pearson correlation coefficient, r ~ -0.38 with a signif¬ 
icance of 0.28, provides no evidence of (linear) correla¬ 
tion between Tkin and N(N2D+)/N(N2H+) ratio. How¬ 
ever, the Spearman rank correlation coefficient, p ~ -0.61 
with a significance of 0.06 (the probability for this sam¬ 
ple to be an uncorrelated system), implies an existence 
of anti-correlation. This features are predicted by chem¬ 
ical models; the N(N2D+)/N(N2H+) ratio is expected to 
increase in cold regions. However, we did not find any 
significant correlation between the N(N2D+)/N(N2H''') 
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Fig. 5.— Line widths of N 2 H+ (1-0) versus N(N2D1’)/N(N2H+) ratios. The colors represent the different type of sources as Fig. 2. The 
error bars indicate Icr uncertainties in our data (red). 

ratio and riHa (Figure 3b). 

To study the evolutionary states of LLOs, we compare 
the chemical evolutionary indicator, N(N2D+)/N(N2H+) 
ratio, and the dynamical evolutionary indicator, 

N2H+ (1-0) line width, with that of star le ss core s 
and Class 0 sourc e s fro m iCraosi et al. I (l2005h ). 
lEmorechtinger et al. I (|2009t) and iFriesen et al. I (|2013 ). 

Figures 4a and 4b show the line width and 
N(N2D+)/N(N2H+) ratio populations, respectively. The 
line widths of Class 0 sources are relatively larger than 
that of starless cores, and most of our LLOs are located 
in the region where Class 0 sources and starless cores 
overlap (Figure 4a). This result implies that LLOs are 
likely young Class 0 sources whereas there is no signif¬ 
icant difference in their N(N2D+)/N(N2H+) ratio (Fig¬ 
ure 5). The N(N2D+)/N(N2H+) ratios (median=0.05) of 
our LLOs are similar to that of starless cores and Class 
0 sources, although only 10 sources have N2D+ (3-2) 
detections. (Though N2D+ J = 3 —>■ 2 is detected in 
DCE 107, the two components are unresolvable, which 
prevents us from deriving the N2D’'' column density for 
each component.) 


5.3. Infall indicators 

Figure 6 shows the comparison of the asymmetry pa¬ 
rameters calculated from HCO+ (3-2) and HCN (3-2). 
The asymmetry parameters obtained from HCO’*' (3-2) 
and HCN (3-2) are quite consistent with each other ex¬ 
cept for DCE 185; HCO+ (3-2) 5v and HCN (3-2) 5v 
show a positive correlation. Eor the 13 targets with both 
HCO+ (3-2) and HCN (3-2) detections, 9 targets have 
the same mathematical sign (Table 6). The remaining 4 
sources are DCE 081, 185, 109 and 063. The asymmetry 
parameters of DCE 109 and 063 have S/N ratios less than 
2(7 for both HCO+ (3-2) and HCN (3-2), which implies 
that no clear systematic motions are occurring (Eigure 
6). The HCO'*' (3-2) spectrum of DCE 185 shows a 
double-peak profile and the peak values are very similar. 


If we used the alternative peak, the asymmetry param¬ 
eter would become -1.12, which would be consistent to 
that from HCN (3-2). However, it is still unclear if DCE 
185 has a clear infall or redshifted signature based on the 
low spectral resolution data. The DCE 081’s HCN (3-2) 
5v indicates no clear systematic motion (^ 1.2a) and its 
HCN (3-2) detection is quite weak (S/N ~ 4.2). Since 
its HCO’*' (3-2) spectrum shows clear asymmetry profile 
(see below), we suggest that its asymmetry parameter 
from HCO'*' (3-2) is much more reliable than that from 
HCN (3-2) for DCE 081. 

Using the two-layer model, the infall velocities of 
13 targets are derived using the HCO+ (3-2) spectra. 
We co mpared the derived inf all velocities with previous 
works. iBelloche et al. I (j2002[) have derived an infall ve¬ 
locity field in DCE 001 (IRAM 04191). The infall veloc¬ 
ity is decreasing from the center and becomes uniform 
(Un 0.10 ± 0.05 km s-i) between -2000-3000 AU 
and —10000-12000 AU. Our derived infall velocity, 0.07 
± 0.03 km s“^ , is consistent with this infall velocity 
field in the outer region. For DCE 004 (L1521E), our 
derived infall velocity 0.14 ± 0.02 is slightly smaller than 
the i nfall velocity 0.2-0.3 km s“^ (2000 to 3000 AU) 
from lOnishi et al. I (|1999f ). but much larger than that 
from iLee et al. I (j2004f l. 0.014 km s“^ and 0. 045 km s~^ 
from C S (2-1) and (3-2) spectra, respectively. iLee et al~l 
(|2004f) concluded that the CS (2-1) and (3-2) may trace 
different regions, which results in a difference between 
the derived infall velocities. HCO+ (3-2) may trace a 
more excited region since it has a higher energy level. 
In addition, the difference of optical depth and deple¬ 
tion between molecules may also explain the discrepancy 
between infall velocity measurements. 

We also compared the infall velocities with the free- 
fall velocities. The free-fall velocities, vs = [GM/ 
l|Keto et al. Il200'il where M and R are the envelope mass 
and radius, respectively, are derived from the RADEX 
best-fit models assuming the central star has a negligible 
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mass. All sources have infall velocities less than free-fall 
velocities except for DCE 038, 063 and 090. The in¬ 
fall velocities of these three sources are not signihcantly 
larger than their free-free fall velocities (within 3 (t un¬ 
certainty of Vin)- Thus, none of our targets has infall 
velocity significantly larger than its free-fall velocity. 

All the sources with positive Tdn have negative HCO+ 
(3-2) 5v and vice versa except for DCE 064. This re¬ 
sult suggests that the infall or red shifted (outward) sig¬ 
natures are quite reliable in the HCO+ infall analysis. 
For DCE 064, the asymmetry spectrum with blueshifted 
wing suggesting an expanding motion is in contrast to the 
negative Sv derived from both HCO+ (3-2) and HCN (3- 
2). The spectral resolution of our current data is insuffi¬ 
cient to identify if DCE 064 has an infalling or expanding 
envelope. We find that the HCO+ spectrum of DCE 064 
has a double-peak profile, and the peak values are similar 
as DCE 185. The HCO"'' dv will be ~ 0.8 (expanding) if 
we use the second peak whereas the 6v derived from the 
first peak is -0.78. In addition, the asymmetry parame¬ 
ter of HCN (3-2), Sv = -0.37, is not significant (S/N ^ 
1.8). Therefore, the blue shifted wing hints that outward 
motion may occur in DCE 064. 

It is noteworthy that the HCO+ systemic velocities 
from the two-layer model fitting are occasionally shifted 
from the systemic velocities obtained from N2H+ (Fig¬ 
ure 2). These results suggest that HCO+ and N2H+ 
trace different regions in the cores, and hint that the 
derived asymmetry parameters could be false. Observa¬ 
tions of optically thin isotopologue H^^CO’*' or HC^®0+ 
are required to correct the infall measurements by pro¬ 
viding an unbiased Vlsr- HCO+ can also trace the 
protostellar outflows which could contaminate our de¬ 
terminations of infall motions. Besides, more compli¬ 
cated dynamic structures such as molecular outflows 
or core rotat i on co uld affect our results. In addition, 
iSmith et al. I (|2012f) found that the blue asymmetry pa¬ 
rameters could be frequently undetected toward collaps¬ 
ing cores l ocated in filaments. Based on radiative trans¬ 
fer models. [Smith et al. I (j20I2ll found that geometry pa¬ 
rameters such as inclination angle will significantly affect 
the profile of optically thick tracers like CS (2-1) and 
HCN (I-O). In ord er to find a reliable infall indicator, 
iChira et al. I (|20I4ll modeled the line profiles of HCO+ 
and HCN among transitions from J = I—J-0toJ = 5 
—>■ 4. They found that J = 4 —)• 3 has the best chance 
to probe infall motions, but the J = 3 —>■ 2 and J = 5 
—>■ 4 transitions of both HCN and HCO’*' are also good 
indicators. 

As a result, 8 sources have positive infall velocities (in¬ 
fall motions) and 5 sources have negative infall velocities 
(outward motions). Out of the 8 infall candidates, 3 
sources (DCE 004, 001, and 092) have significant blue 
asymmetry (S/N > 3) from both HCO’^' (3-2) and HCN 
(3-2). These 3 sources have very low kinetic tempera¬ 
tures (6.5 K < Tkin < 8.0 K) and bolometric tempera¬ 
tures (20 K < Tboi < 27 K). Four sources (DCE 065, 038, 
182, and 090) have a blue asymmetry parameter in ei¬ 
ther HCO’*' (3-2) or HCN (3-2) in the range 2-3 ct, which 
suggests that infall motions may be occurring in these 
candidates too. DCE 065 and 038 also have low kinetic 
temperatures (7.5 and 9.0 K) and bolometric tempera¬ 
tures (29 and 66 K). DCE 182 and 090 have relatively 
high bolometric temperatures (105 and II4K, i.e. young 



Sv HCO+(3-2) 

Fig. 6.— Asymmetry parameters derived from HCO”*" (3-2) 
versus that from HCN (3—2). The blue and red colors indicate the 
area for infall motion and outward motion, respectively. The error 
bars indicate Icr uncertainties. 

Class I sources), and DCE 090 has a high kinetic tem¬ 
perature of 30.5 K (Tkin of DCE 182 cannot be solved 
due to lack of N2H+ (3-2)). Our results suggest that the 
infall motions are likely to occur at an early evolutionary 
stage but can still appear at a later stage ( Class I), which 
is consistent with iGregersen et al. I (|2000l) . 

6. CONCLUSION 

We observed N2H+ (I-O), N2H+ (3-2), N2D+ (3-2), 
HCO"*" (3-2) and HCN (3-2) toward 16 LLOs using the 
KP 12m telescope and SMT in ARO. We derived the 
kinetic temperatures, N2H“'' column densities, H2 densi¬ 
ties and core sizes by fitting spectra with the non-LTE 
code RADEX. We obtain these physical parameters for 
13 LLOs from N2H+ (I-O) and (3-2) observations. We 
derived the N(N2D''')/N(N2H+) ratio which traces the 
evolutionary state of starless cores and protostellar cores. 
To further probe the infall motions, we used HCO+ (3- 
2) and HCN (3-2) as optically thick tracers and N2H“'" 
(I-O) as optically thin tracer. Our conclusions are the 
following: 

1. LLOs tend to have relatively low H2 den- 
sities compared with the Class 0 objects in 
lEmorechtinger et al. I (120(191) . Only DCE 001 and 
185 have higher H2 densities than Class 0 objects; 
however, their densities are most likely overesti¬ 
mated due to the beam size difference in N2H+ 
(I-O) and (3-2) observations. 

2. For Tkin < 10 K, LLOs have N2D+/N2H+ col¬ 
umn density ratios between 0 .03-0.17 similar to the 
Class 0 ob jects observed by lErnorechtinger et al. I 
(|2009l) and lFriesen et al. 1 (I20I3D . and for Tvi„ > 10 
K, LLOs have consistently low N2D+/N2H+ col¬ 
umn density ratios ^ 0.04. The Spearman rank 
correlation test implies an anti-correlation between 
the N2D+/N2H+ column densities and the kinetic 
temperatures in LLOs with a significance of 0.06. 

3. We find that the N2D+/N2H+ column density ra¬ 
tios of LLOs are similar to that of starless cores 
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and Class 0 objects, but the line widths of LLOs 
are mostly populated in the starless cores and Class 
0 objects overlapping region. We suggest that the 
molecular line width could be a better age indi¬ 
cator than N(N2D+)/N(N2H+) ratio at the early 
stage and LLOs are likely young Class 0 protostel- 
lar sources. 

4. We identify eight infall candidates by fitting the 
HCO'*' (3-2) spectra with a two-layer model. Seven 
sources are supported by their asymmetry parame¬ 
ters with S/N > 2 in either or both HCO’*' (3-2) or 
HCN (3-2). Our result suggests that infall motions 
tend to occur at a very early evolutionary stage but 


can still occur in young Class I sources (DCE 090 
and 182). 
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Fig. A1.— The slice in the model grid for T^in versus The black line indicates the confidence interval which fits to the synthetic 

spectra. The color area indicates the models fit to the observing parameters, Tex and r, within their errors. 
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APPENDIX 

We evaluated the feasibility of using N 2 H+ (1-0) and (3-2) data to well constrain the physical parameters with 
RADEX. The method is to find the relation between the uncertainties of Tkin and nuj and a set of synthetic spectra 
with different rms noise levels. We first construct a reference model using RADEX and compute the synthetic N 2 H+ 
(1-0) and (3-2) spectra for this reference model. Then we added artificial noise to the spectra and derived the 
excitation temperatures and optical depths using the HFS fitting method. The data with lower rms noise levels would 
have lower uncertainties in their excitation temperatures and optical depths. We made grid of models with varied 
kinetic temperatures, H 2 densities, and N 2 H+ column densities. Then we compute the for tho reference model and 
each model in the grid based on the excitation temperatures and optical depths. 

Here we describe the process and results in detail. Using RADEX without hyperfine structure, we derive Tex and r 
for both N 2 H+ (1-0) and (3-2) from a reference model set, Tkin = 8.75 K, log(nH 2 ) = 5.225, N(N 2 H+) = 1.1 x 10^^ 
cm“^ and core size = 75". We constructed the synthetic N 2 H+ (1-0) and (3-2) spectra with 

7MB(r') = 4'(J(Tex) - J(Tbg))(l - exp{-T{v))), (1) 

considering the beam sizes of 67.5" for N 2 H+ (1-0) and 27" for N 2 H+ (3-2). We added noise into the synthetic 
spectra. The noise levels are from 0.04 to 0.16 K (at a spectral resolution of 0.078 km s“^) with a step of 0.02 K 
for N 2 H+ (1-0) and from 0.01 to 0.07 K (at a spectral resolution 0.27 km s“^) with a step of 0.01 K for N 2 H+ 
(3-2); in total, seven sets of synthetic spectra are constructed. (The rms noise level of our observational data are 
^0.1 K and ^0.04 K for these two transitions.) We obtain Tex and r using HFS fits to the synthetic spectra. 
As expected, the errors of Tex and r are correlated to the artificial rms noise levels. We then construct a grid of 

models with varied kinetic temperatures, H 2 densities, and N 2 H+ column densities. We calculate the using 

the Tex and r from the synthetic spectra of the reference model and that from each model in the grid, i.e. Ax^ = 
/ rex,l-0-r,x,l-0,model ^2 _|_ / Ti -p -Ti_o,,nodei N 2 / Tex,3-2-T,x.3-2,model \2 _|_ / T 3 - 2 -T3_2,model \2 Al ShOWS the 68.3% 

‘^T'ex,l-0 ‘’’Tex,3-2 ‘^’'3-2 

confidence interval in the slices of Tkin vs. nH 2 space. The color regions represent the models fitting the observing 
parameters within the errors. The regions become broader when the uncertainties become larger with larger rms noise 
levels. With a fixed Tex or r, the fitting models are always populated from high Tkin with low nH 2 to low Tkin with 
high nH 2 - This phenomenon explains that the models with high Tkin low and low Tkin high nH 2 are sometimes 
discriminated with difficulty. As a result, we provided the reference values of noise for deriving the physical parameters. 
For example, the data set with best sensitivity, cfi-q ^ 0.04 K and cr 3_2 0.01 K, can be used to constrain the T^n 

in a range of ^ 7-11 K and the nH 2 in a range of ~ 6.3 x 10^ - 1.6 x 10® cm“® (Figure Al). Note that if we fit the 
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Fig. A2.— Comparison of HFS fitting (blue) and RADEX fitting (red) for DCE 081. The left and right panels show the N 2 H+ (1-0) 
and (3-2) spectra of DCE 081, respectively. The bottom panels show the residuals of the RADEX fit. 

spectra with HSF directly, the degeneracy can be marginally reduced (see §4.1). Our analysis provides the reference 
sensitivity which is required for constraining the physical conditions with RADEX. However, the analysis does not 
guarantee that the sensitivity allows to achieve the constraint as our example. If the best-fit model is located in a 
different region, the confidence interval is expected to be different. 




















